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Abstract 

Various applications of quantum algebraic techniques in nuclear structure physics and 
in molecular physics are briefly reviewed and a recent application of these techniques to the 
structure of atomic clusters is discussed in more detail. 

1. Introduction 

Quantum algebras (also called quantum groups) [|l], H; il are deformed versions of the 
usual Lie algebras, to which they reduce when the deformation parameter q is set equal 
to unity. From the mathematical point of view they are Hopf algebras. Their use in 
physics became popular with the introduction ||^, ^ of the g-deformed harmonic oscillator 
as a tool for providing a boson realization of the quantum algebra sUg(2), although similar 
mathematical structures had already been known Initially used for solving the quantum 
Yang-Baxter equation, quantum algebras have subsequently found applications in several 
branches of physics, as, for example, in the description of spin chains, squeezed states 
10], hydrogen atom and hydrogen-like spectra p, [10|, rotational and vibrational nuclear 
and molecular spectra and in conformal field theories. By now much work has been done 
ri| , |l^ on the g-deformed oscillator and its relativistic extensions [jl5|, |16[, and 

several kinds of generalized deformed oscillators fl^, [l^, ^ and generalized deformed 



su(2) algebras |2^, have been introduced. 

Here we shall confine ourselves to applications of quantum algebras in nuclear structure 
physics and in molecular physics. The purpose of this short review is to provide the reader 
with references for further reading. In addition a recent application of quantum algebraic 
techniques to the structure of atomic clusters will be discussed in more detail. 

2. The suy(2) rotator model 

The first application of quantum algebras in nuclear physics was the use of the de- 



formed algebra sUg(2) for the description of the rotational spectra of deformed and 
superdeformed nuclei. The Hamiltonian of the g-deformed rotator is proportional to 
the second order Casimir operator of the sUq(2) algebra. Its Taylor expansion contains 
powers of J( J+ 1) (where J is the angular momentum), being similar to the expansion 
provided by the Variable Moment of Inertia (VMI) model. Furthermore, the deformation 
parameter r (with q = e*"^) has been found p6| to correspond to the softness parameter 
of the VMI model. Through a comparison of the sUq(2) model to the hybrid model the 
deformation parameter r has also been connected to the number of valence nucleon pairs 



28[] and to the nuclear deformation (3 |2^. Since r is an indicator of deviation from the 
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pure su(2) symmetry, it is not surprising that r decreases with increasing /? [^. The sug(2) 



model has been recently extended to excited (beta and gamma) bands 

B(E2) transition probabilities have also been described in this framework [^. In this 
case the g-deformed Clebsch-Gordan coefficients are used instead of the normal ones. (It 
should be noticed that the g-deformed angular momentum theory has already been much de- 
veloped 1^.) The model predicts an increase of the B(E2) values with angular momentum, 
while the rigid rotator model predicts saturation. Some experimental results supporting 



this prediction already exist |]^. Similarly increasing B(E2) values are predicted by a 
modified version |^ of the su(3) limit of the Interacting Boson Model (IBM), by the su(3) 



limit of the sdg Interacting Boson Model |^ , by the Fermion Dynamical Symmetry Model 



(FDSM) as well as by the recent systematics of Zamffi and Casten ||35|| . 
3. Extensions of the suq(2) model 

The sUg(2) model has been successful in describing rotational nuclear spectra. For 



the description of vibrational and transitional nuclear spectra it has been found |^ that 
J {J + 1) has to be replaced by J {J + c). The additional parameter c allows for the 
description of nuclear anharmonicities in a way similar to that of the Interacting Boson 
Model (IBM) and the Generalized Variable Moment of Inertia (GVMI) model [gV]. The 
use of J( J + c) instead of J( J + 1) for vibrational and transitional nuclei is also supported 
by recent systematics [^. 

Another generalization is based on the use of the deformed algebra su$(2) ||21], 
23i [2^ , which is characterized by a structure function $. The usual su(2) and sug(2) 



algebras are obtained for specific choices of the structure function $. The su$(2) algebra 
has been constructed so that its representation theory resembles as much as possible the 
representation theory of the usual su(2) algebra. Using this technique one can construct, 
for example, a rotator having the same spectrum as the one given by the Holmberg-Lipas 
formula A two-parameter generalization of the sUq(2) model, labelled as sUgp(2), has 
also been successfully used for the description of superdeformed nuclear bands |^ . 

4. Pairing correlations 



It has been found that correlated fermion pairs coupled to zero angular momentum 
in a single-j shell behave approximately as suitably defined g-deformed bosons. After 
performing the same boson mapping to a simple pairing Hamiltonian, one sees that the 
pairing energies are also correctly reproduced up to the same order. The deformation 
parameter used (r = In q) is found to be inversely proportional to the size of the shell, thus 
serving as a small parameter. 

The above mentioned system of correlated fermion pairs can be described exactly by 



suitably defined generalized deformed bosons Then both the commutation relations 



are satisfied exactly and the pairing energies are reproduced exactly. The spectrum of 
the appropriate generalized deformed oscillator corresponds, up to ffist order perturbation 
theory, to a harmonic oscillator with an perturbation. 

If one considers, in addition to the pairs coupled to zero angular momentum, pairs cou- 
pled to non-zero angular momenta, one finds that an approximate description in terms of 
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two suitably defined g-oscillators (one describing the J = pairs and the other correspond- 
ing to the J 7^ pairs) occurs The additional terms introduced by the deformation 
have been found |^ to improve the description of the neutron pair separation energies of 
the Sn isotopes, with no extra parameter introduced. 

g-deformed versions of the pairing theory have also been given in [44, 4^ . 



5. g-deformed versions of nuclear models 

A g-deformed version of a two dimensional toy Interacting Boson Model (IBM) with 
sUg(3) overall symmetry has been developed mainly for testing the ways in which 



spectra and transition probabilities are influenced by the g-deformation. The question of 
possible complete breaking of the symmetry through g-deformation, i.e. the transition from 
the sug(2) limiting symmetry to the sOg(3) one has been examined |^8|, It has been 



found that such a transition is possible for complex values of the parameter g [^. (For 
problems arising when using complex g values see |50] ) . Complete breaking of the symmetry 



has also been considered in the framework of an sUg(2) model |^T|. It has also been found 



that g-deformation leads (for specific range of values of the deformation parameter r, 
with g = e*"^) to a recovery of the u(3) symmetry in the framework of a simple Nilsson 
model including a spin-orbit term. Finally, the Og(3) limit of the toy IBM model has been 
used for the description of ^^O + a cluster states in ^°Ne, with positive results . 



g-deformed versions of the o(6) and u(5) limits of the full IBM have been discussed in 
5^ , p5| [56[] . The g-deformation of the su(3) limit of IBM is a formidable problem, since 



the sug(3) D sOg(3) decomposition has for the moment been achieved only for completely 



symmetric sug(3) irreducible representations [0, 0, |59 



Furthermore a g-deformed version of the Moszkowski model has been developed ^ 



and RPA modes have been studied |£2| in it. A g-deformed Moszkowski model with cranking 
has also been studied ||63| in the mean-field approximation. It has been seen that the residual 
interaction simulated by the g-deformation is felt more strongly by states with large J^. The 
possibility of using g-deformation in assimilating temperature effects is receiving attention. 



since it has also been found that this approach can be used in describing thermal effects 
in the framework of a g-deformed Thouless model for supercoductivity. 

In addition, g-deformed versions of the Lipkin-Meshkov-Glick (LMG) model have been 
developed, both for the 2-level version of the model in terms of an sug(2) algebra |^ 
for the 3-level version of the model in terms of an sUg(3) algebra |]66| . 



and 



6. Anisotropic quantum harmonic oscillator with rational ratios of frequen- 
cies 

The symmetries of the 3-dimensional anisotropic quantum harmonic oscillator with 
rational ratios of frequencies (RHO) are of high current interest in nuclear physics, since they 
are the basic symmetries underlying the structure of superdeformed and hyperdeformed 
nuclei p7| , |6^] . The 2-dimensional RHO is also of interest, in connection with "pancake" 
i.e. very oblate nuclei. Cluster configurations in light nuclei can also be described 
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in terms of RHO symmetries, which underlie the geometrical structure of the Bloch- 



Brink a-cluster model |j7T| . The 3-dim RHO is also of interest for the interpretation of the 
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observed shell structure in atomic clusters ||72|, [73[], especially after the realization that large 
deformations can occur in such systems |7^. (See section 9 for further discussion of atomic 
clusters.) 

The two-dimensional and three-dimensional |]75[ anisotropic harmonic oscillators have 
been the subject of several investigations, both at the classical and the quantum mechanical 
level (see for references). These oscillators are examples of superintegrable systems. 

The special cases with frequency ratios 1:2 and 1:3 have also been considered ||78|. While 



at the classical level it is clear that the su(N) or sp(2N,R) algebras can be used for the 
description of the N-dimensional anisotropic oscillator, the situation at the quantum level, 
even in the two-dimensional case, is not as simple. It has been proved that a generalized 
deformed u(2) algebra is the symmetry algebra of the two-dimensional anisotropic quantum 
harmonic oscillator which is the oscillator describing the single-particle level spectrum 
of "pancake" nuclei, i.e. of triaxially deformed nuclei with » tOy, u^. Furthermore, 
a generalized deformed u(3) algebra turns out to be the symmetry algebra of the three- 
dimensional RHO []77| . 

7. Three-dimensional g-deformed (isotropic) harmonic oscillator 

Recently the 3-dimensional g-deformed (isotropic) harmonic oscillator has been studied 



in detail 



58, 5S . It turns out that 



following the mathematical developments of | 
in this framework, one can reproduce level schemes similar to the ones occuring in the 
modified harmonic oscillator model, first suggested by Nilsson |^0|, An appropriate 
g-deformed spin-orbit interaction term has also been developed Including this term 
in the 3-dimensional g-deformed (isotropic) harmonic oscillator scheme one can reproduce 
level schemes similar to these provided by the modified harmonic oscillator with spin-orbit 
interaction. It is expected that this scheme, without the spin-orbit interaction term, will 
be appropriate for describing the magic numbers occuring in the various kinds of atomic 
clusters |7^, |73|, since a description of magic numbers of atomic clusters in terms of a Nilsson 
model without a spin-orbit interaction has already been attempted |8^. This subject will 
be discussed in some detail in Section 9. 

A recent review of the applications of quantum algebraic techniques to nuclear structure 
problems can be found in |]83| . 

8. The use of quantum algebras in molecular structure 

Similar techniques can be applied in describing properties of diatomic and polytomic 
molecules. A brief list will be given here. 

1) Rotational spectra of diatomic molecules have been described in terms of the sUg(2) 
model [0. As in the case of nuclei, g is a phase factor [q = e*"^). In molecules r is 
of the order of 0.01. The use of the sug(2) symmetry leads to a partial summation of 
the Dunham expansion describing the rotational-vibrational spectra of diatomic molecules 
8^ . Molecular backbending (bandcrossing) has also been described in this framework . 
Rotational spectra of symmetric top molecules have also been considered |]86| , |87[| in the 



framework of the sUq(2) symmetry. Furthermore, two g-deformed rotators with slightly 
different parameter values have been used for the description of A/ = 1 staggering 
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effects in rotational bands of diatomic molecules. (For a discussion of A/ = 2 staggering 
effects in diatomic molecules see |p9| ). 

2) Vibrational spectra of diatomic molecules have been described in terms of g-deformed 
anharmonic oscillators having the sUg(l,l) ||90| or the Ug(2) D Og{2) |]91| symmetry, as well as 



in terms of generalized deformed oscillators similar to the ones described in sec. 3 [92, p3 
These results, combined with 1), lead to the full summation of the Dunham expansion 

qp{2) D Ogp(2) symmetry 



90, 91 



A two-parameter deformed anharmonic oscillator with Ug 



has also been considered ^ ^ 

3) The physical content of the anharmonic oscillators mentioned in 2) has been clarified 
by constructing WKB equivalent potentials (WKB-EPs) and classical equivalent potentials 
|p5[| providing approximately the same spectrum. The results have been corroborated by the 
study of the relation between sUg(l,l) and the anharmonic oscillator with anharminicities 



p6[] . Furthermore the WKB-EP corresponding to the sUg(l,l) anharmonic oscillator has 
been connected to a class of Quasi- Exactly Soluble Potentials (QESPs) . 

4) Generalized deformed oscillators giving the same spectrum as the Morse potential |^ 
and the modified Poschl-Teller potential |9^, as well as a deformed oscillator containing 
them as special cases ||10U | have also been constructed. In addition, (/-deformed versions 



of the Morse potential have been given, either by using the sOg(2,l) symmetry ||101|| or by 
solving a g-deformed Schrodinger equation for the usual Morse potential ||102|| . For the sake 
of completeness it should be mentioned that a deformed oscillator giving the same spectrum 
as the Coulomb potential has also been constructed [|103 |. 

5) A g-deformed version of the vibron model for diatomic molecules has been constructed 
104|| , in a way similar to that described in sec. 5. 



6) For vibrational spectra of polyatomic molecules a model of n coupled generalized 
deformed oscillators has been built |p.05|] , containing the approach of lachello and Oss ||106 



as a special case. In addition a model of two Q-deformed oscillators coupled so that the 
total Hamiltonian has the suq(2) symmetry has been proved | 107 | to be equivalent, to 



lowest order approximation, to a system of two identical Morse oscillators coupled by the 
cross-anharmonicity usually used empirically in describing vibrational spectra of diatomic 
molecules. 

7) Quasi-molecular resonances in the systems ^^C+^^C and have been described 

in terms of a g-deformed oscillator plus a rigid rotator [lOS . 

A review of several of the above topics, concerning the applications of quantum al- 
gebraic techniques to molecular structure, accompanied by a detailed and self-contained 
introduction to quantum algebras, has been given by Raychev ||109|. 



9. The 3-dimensional g-deformed harmonic oscillator and magic numbers of 
alkali metal clusters 



Metal clusters have been recently the subject of many investigations (see [^, |73|, |110 
for relevant reviews). One of the first fascinating findings in their study was the appearance 
of magic numbers ||1 1 1| , |112| , p.l3| , |114| , |115| , |116| , |117|| , analogous to but different from the 
magic numbers appearing in the shell structure of atomic nuclei ||118|| . This analogy led 



to the early description of metal clusters in terms of the Nilsson-Clemenger model [R^ 
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which is a simphfied version of the Nilsson model [BD], PT| of atomic nuclei, in which no 
spin-orbit interaction is included. Further theoretical investigations in terms of the jellium 
model ||119| , |12CI|| demonstrated that the mean field potential in the case of simple metal 
clusters bears great similarities to the Woods-Saxon potential of atomic nuclei, with a slight 
modification of the "wine bottle" type | 121 |. The Woods-Saxon potential itself looks like 
a harmonic oscillator truncated at a certain energy value and flattened at the bottom. It 
should also be recalled that an early schematic explanation of the magic numbers of metallic 
clusters has been given in terms of a scheme intermediate between the level scheme of the 3- 
dimensional harmonic oscillator and the square well . Again in this case the intermediate 



potential resembles a harmonic oscillator flattened at the bottom. 

On the other hand, modified versions of harmonic oscillators Q have been recently in- 
vestigated, as it has already been mentioned. The spectra of g-deformed oscillators increase 
either less rapidly (for q being a phase factor, i.e. q = e*"^ with r being real) or more rapidly 
(for q being real, i.e. q = with r being real) in comparison to the equidistant spectrum 
of the usual harmonic oscillator |53|, while the corresponding (WKB-equivalent) potentials 
p5[| resemble the harmonic oscillator potential, truncated at a certain energy (for q being 



a phase factor) or not (for q being real), the deformation inflicting an overall widening or 
narrowing of the potential, depending on the value of the deformation parameter q. 

Very recently, a g-deformed version of the 3-dimensional harmonic oscillator has been 
constructed |7^, taking advantage of the Ug(3) D sOg(3) symmetry ||122| , p7|] . As it has al- 



ready been mentioned in Section 7, the spectrum of this 3-dimensional g-deformed harmonic 
oscillator has been found [79| to reproduce very well the spectrum of the modified harmonic 
oscillator introduced by Nilsson [^, |8ll|, without the spin-orbit interaction term. Since the 
Nilsson model without the spin orbit term is essentially the Nilsson-Clemenger model used 
for the description of metallic clusters [Q, it is worth examining if the 3-dimensional q- 
deformed harmonic oscillator can reproduce the magic numbers of simple metallic clusters. 
This is the subject of the present section. 

The space of the 3-dimensional g-deformed harmonic oscillator consists of the com- 
pletely symmetric irreducible representations of the quantum algebra Ug(3). In this space 
a deformed angular momentum algebra, sOg(3), can be defined The Hamiltonian of 



the 3-dimensional g-deformed harmonic oscillator is defined so that it satisfies the following 
requirements: 

a) It is an sOg(3) scalar, i.e. the energy is simultaneously measurable with the g-deformed 
angular momentum related to the algebra sOg(3) and its ^-projection. 

b) It conserves the number of bosons, in terms of which the quantum algebras Ug(3) and 
sOg(3) are realized. 

c) In the limit g ^ 1 it is in agreement with the Hamiltonian of the usual 3-dimensional 
harmonic oscillator. 

It has been proved ||79|] that the Hamiltonian of the 3-dimensional g-deformed harmonic 
oscillator satisfying the above requirements takes the form 



H, = huo [N]q 



N+l 



[2] 



where is the number operator and C^^^ is the second order Casimir operator of the algebra 
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sOg(3), while 

M = ^ (2) 

is the definition of g-numbers and g-operators. 

The energy eigenvalues of the 3-dimensional g-deformed harmonic oscillator are then 



79 



E,{n, I) = huo |Ng"+i - + l]| , (3) 

where n is the number of vibrational quanta and / is the eigenvalue of the angular momen- 
tum, obtaining the values I = n,n — 2, ... ,0 or 1. 

In the limit of g — > 1 one obtains limq^iEq{n,l) = huJon, which coincides with the 
classical result. 

For small values of the deformation parameter r (where g = e"^) one can expand eq. (3) 
in powers of r obtaining 



Eq{n, I) = hujQin — huoT{l{l + 1) — n{n + 1)) 

- hiUoT'' + 1) - in(n + l)(2n + 1)) + Oir^). (4) 
The last expression to leading order bears great similarity to the modified harmonic 



oscillator suggested by Nilsson ^ (with the spin-orbit term omitted) 



\ f L 



where 



2. iV(iV + 3) 



< V >^= (6) 

It has been proved that the spectrum of the 3-dimensional g-deformed harmonic os- 
cillator closely reproduces the spectrum of the modified harmonic oscillator of Nilsson. In 
both cases the effect of the l{l + 1) term is to flatten the bottom of the harmonic oscillator 
potential, thus making it to resemble the Woods-Saxon potential. 

The level scheme of the 3-dimensional g-deformed harmonic oscillator (for huQ = 1 and 
T = 0.038) is given in Table 1 ||123|| , up to a certain energy. Each level is characterized 
by the quantum numbers n (number of vibrational quanta) and / (angular momentum). 
Next to each level its energy, the number of particles it can accommodate (which is equal 
to 2(2/ + 1)) and the total number of particles up to and including this level are given. If 
the energy difference between two successive levels is larger than 0.39, it is considered as a 
gap separating two successive shells and the energy difference is reported between the two 
levels. In this way magic numbers can be easily read in the table: they are the numbers 
appearing above the gaps, written in boldface characters. 

The magic numbers provided by the 3-dimensional g-deformed harmonic oscillator in 



Table 1 are compared to available experimental data for Na clusters ||1 1 1| , |112| , |113| , |115| , |116 
in Table 2 (columns 2-6) ||123|| . The following comments apply: 



7 



i) Only magic numbers up to 1500 are reported, since it is known that filling of electronic 
shells is expected to occur only up to this limit ||111|| . For large clusters beyond this point 



it is known that magic numbers can be explained by the completion of icosahedral or 
cuboctahedral shells of atoms 



111 



ii) Up to 600 particles there is consistency among the various experiments and between 
the experimental results in one hand and our findings in the other. 

iii) Beyond 600 particles the predictions of the three experiments, which report magic 
numbers in this region, are quite different. However, the predictions of all three experiments 
are well accommodated by the present model. Magic numbers 694, 832, 1012 are supported 
by the findings of both Martin et al. 



11 111 and Brechignac et al. [p.16 



magic numbers 
magic 



111 



1206, 1410 are in agreement with the experimental findings of Martin et al. 
numbers 912, 1284 are supported by the findings of Brechignac et al., while magic numbers 
676, 1100, 1314, 1502 are in agreement with the experimental findings of Pedersen et al. 



115 |. 



In Table 2 the predictions of three simple theoretical models [ [118| ] (non-deformed 3- 
dimensional harmonic oscillator (column 9), square well potential (column 8), rounded 
square well potential (intermediate between the previous two, column 7) ) are also reported 
for comparison. It is clear that the predictions of the non-deformed 3-dimensional harmonic 
oscillator are in agreement with the experimental data only up to magic number 40, while 
the other two models give correctly a few more magic numbers (58, 92, 138), although they 
already fail by predicting magic numbers at 68, 70, 106, 112, 156, which are not observed. 

It should be noticed at this point that the first few magic numbers of alkali clusters 
(up to 92) can be correctly reproduced by the assumption of the formation of shells of 
atoms instead of shells of delocalized electrons ||124|| , this assumption being applicable 
under conditions not favoring delocalization of the valence electrons of alkali atoms. 

Comparisons among the present results, experimental data (by Martin et al. | |1 1 1| | 
(column 2), Pedersen et al. ||115|| (column 3) and Brechignac et al. ||116|| (column 4) ) 



and theoretical predictions more sophisticated than these reported in Table 2, can be made 
in Table 3 |[L23|| , where magic numbers predicted by various jellium model calculations 



(columns 5-8, ||111| , |112| , |74[), Woods-Saxon and wine bottle potentials (column 9, 



||125|| ), as well as by a classification scheme using the 3n + / pseudo quantum number 
(column 10, ||111|| ) are reported. The following observations can be made: 

i) All magic numbers predicted by the 3-dimensional g-deformed harmonic oscillator are 
supported by at least one experiment, with no exception. 

ii) Some of the jellium models, as well as the 3n + I classification scheme, predict magic 
numbers at 186, 540/542, which are not supported by experiment. Some jellium models 
also predict a magic number at 748 or 758, again without support from experiment. The 
Woods-Saxon and wine bottle potentials of Ref. [|1 25|| predict a magic number at 68, 
for which no experimental support exists. The present scheme avoids problems at these 
numbers. It should be noticed, however, that in the cases of 186 and 542 the energy gap 
following them in the present scheme is 0.329 and 0.325 respectively (see Table 1), i.e. 
quite close to the threshold of 0.39 which we have considered as the minimum energy gap 
separating different shells. One could therefore qualitatively remark that 186 and 542 are 
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"built in" the present scheme as "secondary" (not very pronounced) magic numbers. 
The following general remarks can also be made: 

i) It is quite remarkable that the 3-dimensional g-deformed harmonic oscillator repro- 
duces the magic numbers at least as accurately as other, more sophisticated, models by 
using only one free parameter {q = e"^). Once the parameter is fixed, the whole spectrum is 
fixed and no further manipulations can be made. This can be considered as evidence that 
the 3-dimensional g-deformed harmonic oscillator owns a symmetry (the Uq(3) D sOg(3) 
symmetry) appropriate for the description of the physical systems under study. 

ii) It has been remarked ||111|| that if n is the number of nodes in the solution of the 
radial Schrodinger equation and / is the angular momentum quantum number, then the 
degeneracy of energy levels of the hydrogen atom characterized by the same r?, + / is due to 
the so (4) symmetry of this system, while the degeneracy of energy levels of the spherical 
harmonic oscillator (i.e. of the 3-dimensional isotropic harmonic oscillator) characterized 
by the same 2n + / is due to the su(3) symmetry of this system. 3n + / has been used 
|1 1 1|| to approximate the magic numbers of alkali metal clusters with some success, but no 
relevant Lie symmetry could be determined. In view of the present findings the lack of 
Lie symmetry related to 3?7, + / is quite clear: the symmetry of the system appears to be a 
quantum algebraic symmetry (ug(3)), which is a nonlinear extension of the Lie symmetry 
u(3). 

iii) An interesting problem is to determine a WKB-equivalent potential giving (within 
this approximation) the same spectrum as the 3-dimensional g-deformed harmonic oscilla- 



tor, using methods similar to these of Ref. [95|. The similarity between the results of the 



present model and these provided by the Woods-Saxon potential (column 9 in Table 3) 
suggests that the answer should be a harmonic oscillator potential flattened at the bottom, 
similar to the Woods-Saxon potential. If such a WKB-equivalent potential will show any 
similarity to a wine bottle shape, as several potentials used for the description of metal 



clusters do [|119| , |12CI| , |121|| , remains to be seen. 

In summary, we have shown in this section that the 3-dimensional g-deformed harmonic 
oscillator with Ug(3) D sOg(3) symmetry correctly predicts all experimentally observed magic 
numbers of alkali metal clusters up to 1500, which is the expected limit of validity for the- 
ories based on the filling of electronic shells. This indicates that Ug(3), which is a nonlinear 
deformation of the u(3) symmetry of the spherical (3-dimensional isotropic) harmonic os- 
cillator, is a good candidate for being the symmetry of systems of alkali metal clusters. 

This work has been supported by the Greek Secretariat of Research and Technology 
under contract PENED 95/1981. 
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Tabic 1: Energy spectrum, Eq{n,l), of the 3-dimensional g-deformed harmonic oscillator 
(eq. (3)), for Huq = 1 and q = with r = 0.038. Each level is characterized by n (the 
number of vibrational quanta) and / (the angular momentum). 2(2/ + 1) represents the 
number of particles each level can accommodate, while under "total" the total number 
of particles up to and including this level is given. Magic numbers, reported in boldface, 
correspond to energy gaps larger than 0.39, reported between the relevant couples of energy 
levels. 
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Table 2: Magic numbers provided by the 3-dimensional g-deformed harmonic oscillator 
(Table 1), reported in column 1, are compared to the experimental data of Martin et al. 
[Ill] (column 2), Bj0rnholm et al. [112] (column 3), Knight et al. [113] (column 4), Pedersen 
et al. [115] (column 5) and Brechignac et al. [116] (column 6), concerning Na clusters. The 
magic numbers provided [118] by the (non-deformed) 3-dimensional harmonic oscillator 
(column 9), the square well potential (column 8) and a rounded square well potential 
intermediate between the previous two (column 7) are also shown for comparison. See text 
for discussion. 
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Table 3: Magic numbers provided by the 3-dimensional g-deformed harmonic oscillator 
(Table 1), reported in colunmn 1, are compared to the experimental data of Martin et al. 
[Ill] (column 2), Pedersen et al. [115] (column 3), and Brechignac et al. [116] (column 
4), as well as to the theoretical predictions of various jellium model calculations reported 
by Martin et al. [Ill] (column 5), Bj0rnholm et al. [112] (column 6), Brack [73] (column 
7), Bulgac and Lewenkopf [74] (column 8), the theoretical predictions of Woods-Saxon and 
wine bottle potentials reported by Nishioka et al. [125] (column 9), as well as to the magic 
numbers predicted by the classification scheme using the 3n + / pseudo quantum number, 
reported by Martin et al. [Ill] (column 10). See text for discussion. 
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